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Abstract 



We study the parabolic free boundary problem of obstacle type 
a du f 



Under the condition that / = Hv for some function v with bounded 
second order spatial derivatives and bounded first order time derivative, 
we establish the same regularity for the solution u. Both the regularity 
and the assumptions are optimal. 

Using this result and assuming that / is Dini continuous, we prove 
that the free boundary is, near so called low energy points, a C 1 graph. 

Our result completes the theory for this type of problems for the heat 
operator. 

1 Introduction 

We present a proof of the interior optimal regularity, that is, W^-regularity 
(bounded second order spatial derivatives and bounded first order time deriva- 
tive) for solutions to the so-called parabolic no-sign obstacle problem 



under minimal assumptions on the data. Here B r is the unit ball, Q~ = B r x 
(— r 2 ,0], / = Hv where v £ W^ 1 . Since we are interested in the interior 
regularity, the assumptions on g are not very important but we assume that g 
is bounded and continuous for the sake of definiteness. 

Naturally, we cannot expect that a solution u of ([1]) is in any better regularity 
class than the solution, v, of the heat equation with right hand side / 



Noticing that if / £ L p then fx{u^o} £ L p so for / in the L p -range Calderon- 
Zygmund theory directly implies that u £ W^ 1 - which is as good regularity as 
v. 




l i 



(1) 



Hv = f. 
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It is also easy to see that u is no better than C 11 — regular (bounded second 
order spatial derivatives) in the space directions even in the case when / is 
constant. As a matter of fact, with an appropriate choice of g the function 

u(x,t) = — (max(0, xi)) 2 

is a solution with f(x,t) = 1. Similarly, the example 

u(x,t) = -(t - 1/2)+ 

shows that u(x, t) is no better than C ' 1 regular (bounded first order time 
derivative) in time. This shows that even for constant / we cannot hope for 
better regularity than u G W^ 1 . This leads to the question: What is the 
weakest possible assumptions on / that assures that a solution u(x,t) of |T]) 
satisfies u £ W^ 1 ? 

Since, in general, a solution to (jTJ) is never more regular than a solution to 
Hv = f one might ask: If / is such that the solution v to Hv — f satisfies 
v e W^ 1 will a solution to also satisfy u G W^ 1 ? 

In the main theorem of this paper we answer the question in the affirma- 
tive. This clearly provides an optimal regularity result for the no-sign parabolic 
obstacle problem. Our main theorem is the following. 

Theorem 1. Let u be a solution to (QJ) and assume furthermore that f = Hv 
in Qi where v G W^ > 1 (Q7) and that u £ L 1 (Q7). Then u G W£ (QT> 2 ) an< ^ 

W D2u h^Q- /2) + IKIU( Ql>) < c(IMI ilW r) + WDM^iQ-) + Hh^Q-)), 

(2) 

where C depends on the dimension. 

It might be in order to comment on the assumption that Hv — f where 
v G W^lQi). The assumption might look technical, but to the authors' 
knowledge there is no necessary condition on / that assures that the solution v 
to Hv = f is in Wli 1 . The weakest sufficient condition known to the authors 

<t CO 

is / G C Dini (Qi) (see Section O and for instance [13] or [16]). Naturally, 
the theorem holds, with minor changes in ([2]), under the assumption that / G 
C uini (Qi). But the assumption we make on / is strictly weaker than Dim- 
continuity and it also highlights that up to W^ 1 the regularity of solutions to 
the no-sign parabolic obstacle problem are as regular as the solutions to the 
corresponding heat equation. 

Before we give a brief sketch of the history of the problem and formulate 
our second main theorem we would like to remark that this is a free boundary 
result. As a matter of fact the problem can be formulated as follows. Given 
/ G L p (Qi) and g G C(d p Q±) find a set fi C and a function u G W^{Q^) 
solving 

Hu = fxn inQr, 
u = |Vu| =0 in 
u = g on dpQi . 
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The free boundary is <9 p interior({M = 0}) D Qi — d p {u ^ 0} D Qi ■ In our proof 
we will use that u = outside of spt(ifu) (to be more precise, we use that 
||_D 2 u|| — u t — a.e.). The above theorem is not true, in general, for solutions 
to 

( Hu = fxs in Q l , 
\ u = g on dpQi, 

where £ is an arbitrary set. A simple counterexample for the time independent 
case, with S = {xiX2 > 0}, is given in [2]. In [2] a solution, u, to (j3|) with 
f(x, t) = —1 and £ = {x\X2 > 0} is explicitly calculated and u(x, t) = u(x, 0) ^ 
C M (.Bi(0)). That /(a:,i) < in 0, whereas we usually think of / > in 
obstacle problems is irrelevant since ([3]) is linear, so we can simply change the 
sign of the equation by taking — u instead of u. The moral sense of the example 
is that the solution to the free boundary problem choses a zero level set that 
assures that the solution has slightly better regularity than what one in general 
would expect a solution to the similar problem ([3]) to have. 

Before we continue we would like to sketch an outline of some previous 
research into this problem. The aim of this outline is not to provide a historical 
survey but to situate our result in the current theory of the obstacle problem. 

The papers we discuss below ([4], [8], [H] [7]) are technically very sophisti- 
cated and we have to refer the reader to the original sources for the full details. 
It should be mentioned that we will, rather mischievously, slightly change the 
conceptual framework of the above papers into the BMO framework of this 
paper in our explanations. 

If / e £°°(Qr) then ii; directly follows that D 2 u e BMO(Q~ /2 ), i.e., 

\\D\- (D\) Q - (x0) \\ L2iQ - (x0)) < Cr^ 2 (|M| Ll(Qr) + ll/IU (Qr) ) , (4) 

where (D 2 u)q- , x0 ) denotes the average over the parabolic half cylinder Q~ (X°) 

and D 2 is the second spatial and first time derivative (see the list of notation 
at the end of the introduction). If \(D 2 u)q-^ xq ^ < C for some constant C 

independent of r and X° then by the triangle inequality \\D 2 u\\ l2 ^q- ^ x a^ < 
C r (™+2)/2 which implies that D 2 u G L°°. It is not difficult to see (cf. Lemma 
[12]) that instead of subtracting (D 2 u)q- ^ x0 ^ in (jH) we can use D 2 p u ^ r ^ x o(x,t) 
where p u , r ,x° t) is a parabolic polynomial that is second order homogeneous 
in x and first order homogeneous in t. In particular it is enough to estimate 
||?Vr-,x IIl°°(q-) m order to derive W^ 1 —regularity for u. 

The first regularity results for parabolic obstacle problems where obtained 
under the assumption that u > 0, which implies that p u ,r,x° ^ ~ C in 
since ||w — p u , r .x° \\l 2 {q-(x )) - Cr 2+ (™+ 2 ^ 2 , by our BMO estimate. But if 
Pu,r,x° > — C and is a second order caloric polynomial then it directly follows 
that p u ,T,x a < C in Qj" which implies the optimal regularity by the above. The 
real difficulties therefore occur for no-sign obstacle problems, i.e., when there is 
no assumption on the sign of the solution. 

The first major breakthrough in the regularity theory for parabolic obstacle 
problems without a sign assumption was achieved in [4], where it was proved 
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that if f{x,t) — 1 then the solution to ([T]) is in W^ 1 . The proof is based on a 
monotonicity formula first proved in [5]. The monotonicity formula is applied 
on the positive and negative parts of D e u, the directional derivatives of u, which 
can be shown to be sub-caloric functions in their supports. The monotonicity 
formula gives uniform bounds, in r, of the following averages 

i J J \V(D e ii) + \ 2 G{x,-s)dxds J J \\7{D e ii)-\ 2 G{x,~s)dxds. (5) 

Since, BMO-estimates implies u = p u ^ } x° up to an error that is bounded in 
it follows that the expression in ([5]) is bounded with p u ^ t x° in place of u. 
Using that p u ^ : x° is a parabolic second order polynomial implies, by elementary 
calculations, that p Uyr ^x° is bounded. 

Since the monotonicity formula in [5] is valid only if (D e u)^ are sub-caloric 
functions the results in [3] are only valid in the case when / is constant. However, 
a refined version of the monotonicity formula, valid when H(D e u)^ > — C for 
some constant C, was proved in [8] . This monotonicity formula makes it possible 
to prove W^ 1 — regularity for solutions when / £ C ' 1 . Clearly, if / ^ C 0,1 then 
i?(£> e u) ± > —C is no longer true which makes it difficult to use this method to 
prove regularity for / less regular than Lipschitz. 

The monotonicity formula approach therefore provides optimal regularity 
results for / £ C 0,1 . There is however a substantial and rather unsatisfying 
gap in the regularity theory. If / £ L°°, then classical methods implies that 
u £ W 2 ' 1 for any p < oo. But in order to achieve the W^ 1 — regularity with the 
above mentioned methods, one needs to assume that V/ £ L°°, i.e., a whole 
extra derivative is required. 

There is another approach to the regularity for the no-sign parabolic obstacle 
problem, which is based on an extra assumption on the behaviour of the free 
boundary. It is not difficult to prove that the solution is W^ 1 close to points X° 
where \{u = 0}nQ~(X°)\ > e\Q~(X°)\ for every r > 0. The most sophisticated 
result of this kind is [Tl] in the elliptic case and [7] for the parabolic case. The 
assumptions on the free boundary are, in order to be as week as possible, rather 
technical so we will have to refer the readers to the original papers for the 
details. For our purposes it is enough to remark that even though the methods 
in [14j and [7] are strong enough to prove optimal regularity of the solution - 
they are only able to do so under assumptions on the free boundary which are 
unfortunately not verifiable in general. 

The proof in our paper is based on the method in 1 where we prove similar 
results for the elliptic problem. Let us sketch the proof of Theorem [1] before we 
state the second main theorem of the paper. If u is a solution to ([I} then we 
can write 

u = v + g + second order caloric polynomial, 

where g is a solution to Hg = — fx{u=o}- We will slightly change the nota- 
tion and write the second order caloric polynomial as S(u, r, X°)p u r x o where 
supg- |p u . ri jfo| = 1, that is, S(u, r, X ) controls the norm of the caloric polyno- 
mial. We will choose S(u, r, X )p u ,r,x° by means of a projection operator (see 
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Definition [TU|) that assures that S(u, r, X°)p urX ° closely approximates u. As 
before, the BMO estimates implies that u S W^ 1 if and only if S(u,r,X°) is 
bounded. 

The idea of the proof is to use that on A = {u = 0} we have D 2 u — and 
thus 

= \\D 2 u\\ L2{KnQ - {x0)) > \S{u,r,X°)\\\b 2 Pu ^ x o\\ L2{AnQ - {x0)) (6) 

~\\D 2 9\\L 2 (AnQ-(x°)) ~ 11-^ u llz 2 (AnQ;:(X<>))- 
Now, since v £ W^ 1 and since \\p u ,r,x°\\L^(Q-) = ■*-) ^ follows that 

II^Vr,Xo|| L »(AnQr(XO)) < \/\ A ^Qr(X )\ 

and 

P 2 «llL 2 (AnQ- ( xo)) < yJ\AnQr(X°)\\\D 2 v\\ L ~. 

Moreover, g may be written g = g + h, where h is a caloric function satisfying, 
by BMO estimates 

||^|| L2(Q , (X0)) < Cr("+ 2 )/ 2 (|| M || Ll(Qr) + WfW^Q-J, 

and g has zero boundary data and thus, by parabolic estimates 



P^ll i2(Q -(Xo)< c ll^ll i2(Q -(xo))<cVI An ^( x °)ill/ll^(Qr)- 

If we disregard the caloric function h then ^ can be written as 

\S(u,r,X°)\ < C||/|| LOO(Qr) , 

which is our desired estimate. We may, unfortunately, not disregard h and this 
explains why the paper is around 20 pages and not just a few lines. 

In order to salvage something out of the above calculation we use the parabolic 
estimate 

II^ )2/i IIl=(q- /2 (x )) - (n+2)/2 ||-C ,2/l llL2(Q-( X o)), 

which means that we can use the calculation ([6]) in Q~, 2 (X°). This estimate 
appears, then applied directly on g — g + h, in (fT4"f and the calculations fol- 
lowing (Q3]) explicates how this revised estimate controls the dyadic decay of 
the measure |A n Qa-j|- This is the heart of the paper and carried out, with 
slight variations, in Proposition[T3] The dyadic decay of the measure \ AnQ 2 -j |, 
Proposition [T4J implies that if S(u,r, X°) is large enough then 

|AnQ- /2 pr°)| l|AnQ-(A°)| 



i.e., the function \A satisfies a Morrey space condition at A - at least at the 
scale r. This is utilised, in Proposition [15J to show that g is indeed small as 
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long as S(u,r,X°) is large in comparison to |M| L i(Q--j and ||-D 2 «||_[,cx>(q-\- For 
full details see the main body of the paper. 

Our method improves on the previous research. It is strong enough to prove 
the optimal regularity with minimal assumptions on / and therefore improves 
on the monotonicity formula approach in [4] or [8] that required / G C ' , We 
do not make any assumptions on the solution or its free boundary (such as 
in [T3] or [?])■ Furthermore, whereas previous methods have utilised powerful, 
but rather specialized, monotonicity formulas with limited reach, our method is 
based on standard IP estimates. This makes it likely that our methods can be 
extended to cover other equations, such as equations with variable coefficients, 
higher order equations or fully non-linear equations. 

As soon as W^ 1 — regularity have been established we are in the position to 
apply the powerful free boundary regularity results developed in [4] , [7] and [12] . 
For that we need to assume that / G C m . This assumption is also optimal in 
the sense that if / g" C Dlnl , then there is a time independent solution for which 
the free boundary is not C 1 at the origin, but in fact a spiral point, see [3]. 
We will only sketch the proof to highlight some minor differences, the reader is 
referred to the original papers for the details. 

If u is a solution to ([1]) then we define the free boundary, in the usual way, 
to be 

T{u) = dpinterior ({u = 0}) , (7) 

where d p stands for the parabolic boundary, referring to the points that are not 
parabolically interior, i.e. the points such that 

Q~ n (interior{u = 0}) = 0, 

for any r small enough. In this terminology our second main theorem is the 
following. 

Theorem 2. Let u be a solution to (J^l and assume in addition that f G 
C Dini (Qr) and f (0,0) = 1. 

• // the origin is a low energy point ( as in Definition \19\) then the free 
boundary T(u) (as defined in ^) is, in a neighbourhood of the origin 
(which might depend on the solution u), a (parabolic) C 1 regular graph. 

• There is a modulus of continuity a andro > (both depending on \\u\\ L1 rQ-< 
and || / 1| QDiKi(Q-)) ) such that if 

MP ({x:u(x,-r 2 ) = 0}nB r ) 

" > (8) 

for some r < ro, then T(u) is a C 1 regular graph in Q~j 2 - Here MD stands 
for the minimal diameter. 

It is in order to explain the assumptions in Theorem [2] for the non-expert 
reader. It is well known that the free boundary is not C 1 everywhere (see [15j for 
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an example in the time independent case) and some extra assumption is needed 
to exclude that the origin is a singular point of the free boundary. A natural 
and correct assumption could be that the zero set of u has positive Lebesgue 
density at the origin, which is a slightly stronger assumption that (jHJ. 

There is also a different statement, based on a monotonicity formula, that 
excludes singularities. The powerful monotonicity formula states that a certain 
energy W(r; u, /, (0, 0)) (defined in Section [6j> is almost increasing in r if u is a 
solution to ([!]) and (0,0) € T(u). Furthermore, W(0 + ;u, /, (0,0)) (which is well 
defined due to the almost monotonicity) can only assume the values 15, 15/2 
or 0, and the value carries geometric information of the free boundary at the 
origin. For reasons explained in [3], the value in does not occur if (0, 0) G T(u). 
If the value is 15/2, we say that the origin is a low energy point. In particular, 
by the discreteness of the limiting energies we can conclude that if © or a 
Lebesgue density condition holds then we are at a low energy point. The energy 
condition we impose is therefore weaker than a Lebesgue density condition but 
comparable to ©. 



1.1 Notation and assumptions 

Throughout the paper we use the following notation: 

X = (x,t),X° = (x ,t°),Y=(y,s) 
u t = d t u=^ 

a. = EIU 

Hu = Au — ut 

Q r {x ,to) = B r (x ) x (t - r 2 ,t + r 2 ) 

Q~(x ,t ) = B r (x ) x (t - r 2 ,0] 

Q r = Q r (0,0), Q~ = Q-(0,0), d p Q- = d p Q-(0,0) 

W^(A) 



D 2 u 



D 2 u 
d t u 



\D 2 u\ = \ D '. 



2 J I 2 J 



\d t u\ 



A = {M = o}nQj; 
r = a p ({ u ^o}°)nQr 



A 



|AnQ~| 
(f)r,X° ~ inzixow II 



rn /?t 



- general points in R™ x K + 

- the time derivative 

- the spatial gradient 

- the Laplace operator 

- the heat operator 

- a parabolic cylinder 

- a half cylinder 

- simplified notation 

- the Sobolev space of functions in 
L P (A) with p integrablc 

second derivatives in space and 
first derivatives in time. 

- the parabolic second derivative 

- the matrix norm used 

- the coincidence set 

- the free boundary, the part of d({u ^ 0}°) 
that is not parabolically interior in A 

- the density of A at the level r 

- the mean value of / over Q~(Xq) 

- the mean value of / over S 



Definition 3 (Dini continuity). A function f(x,t) is said to be Dini continuous 
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(in the parabolic setting) if 



\f(x,t) - f(y, s)\ < a(y/\x - y\ 2 + \t- s\), 
where a is a non-negative continuous function such that <r(0) = and 



J 

Jo 



-s < 00. 



2 Parabolic spaces 

Here we present some useful definitions and result for parabolic spaces. We will 
use the standard notation 

for the average of a function / over a set S. 

Definition 4 (Parabolic BMO). We say that a function f G L 2 (fl) is in 
BMO{VL) if 

imi2 _ „„„ 1 f if/., „\ / r\ 12,11^112 



\\J\\bmo(q)= SU P in , rt | / l/(l/»*)-(/)or(x,t)l +ll/lli=(n) < °°- 

The result below is well known and can be found in for instance [5J. 

Theorem 5 (BMO-estimates for the heat equation). Let Hw = f in Q^. If 

f e L°°{Q R ) then 

W D2w \\bMO(Q r/2 ) + \\ d M\ B MO(Q R/2 ) < C (Wfh~(Q R ) + MlHQ-))' 

Here the constant depends only on the space dimension. 

We will need the following standard results: 
Lemma 6. Let u be a solution of Hu = in Q~ . Then 

C 

W^ 2u Wl^(q- /2 ) < ^W u \\LHQ-y 

The proof of this lemma is contained in the proof of Theorem 8 on page 59 
in 0. 

Lemma 7. Let 

Hu = f in Q~ 
u = on d p Q~ 

then 

||^|| L2(Q , ) <Cr("+ 2 )/ 2 ||/|| i2(Q - ) . 
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This is a rescaled version of Corollary 7.16 in 10 . 
Lemma 8. Let 

Hw = f in Qi , 

where f — Hv for some v satisfying \\D 2 v\\La a ^Q 1 j < oo. Then there exists a 
constant C such that 

\\D 2 w\\^ (Q - /2) < C (\\w\\ LHQT) + ||^||^ (Qn ) • 

Below follows a parabolic version of the Poincare inequality. 

Lemma 9 (Parabolic Poincare inequality). Assume that w £ L 1 (Q^ (0)) . Then 
for some n S (7/8, 1) and some C > 

T - H Q - - x ■ {Vw) Q - | l2(q _ ) < C (||I> a «>|| ia(<3r) + WMvuq-)) ■ 

Proof. We notice that from the Poincare inequality applied of each i-section of 
w we may deduce that 



\ w — (w)B K x{t} — x • (Vw)b k x{*}| dxdt < C / |_D 2 ui| dxdt. (9) 

JQ, 

Hence, it suffices to find proper estimates for the differences 

IIHs«x{t} - ( w )q-Wl2(q-)> IK ( Vw )s K x{t} -x- (Vw) Q -\\ L2{Q - y 

For that purpose we first of all remark that from the Poincare inequality em- 
ployed for fixed x one can conclude 



\w - H ( _ K 2, 0) |^xt < c / Krxj. (10) 
From Holder's inequality we can then conclude 

IHs K x{ t } - H Q -| 2 xt < 

<C [ ([ \w-(w) { _ K2 , 0) \ 2 )xt<C [ \w t \ 2 xt. 
Jq- \Jb k J Jq- 

Now we compute, and use for the third equality that x-V y w(y, ■) = div y (xw(y, 
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together with the divergence theorem 



Q,7 



x ■ (Vio) BsX(t } - x ■ (Vio) ( 
1 



xt 



1 



B. 



V y w(y,t)dy - x ■ 



x ■ vw(y, t)dy 



S7 y w{y, s)dyds 



xt 



K \B K 



k\ J-k 2 JOB 



x ■ vw{y, s)dyds 



xt 



<C / \w(y,t) - (w(y,s)) se (- K 2 >0 )\ yxt 
<C / \w(y,t) - (w(y, s)) se( _ K 2. 0) | 2 t. 

J-k, 2 JdB K 

Using polar coordinates and the mean value theorem, there exist a k G (7/8, 1) 
such that 



\w{y,t) - (io(y,s)) ae (_ K a )0 )| yt 



< C / / \w(y,t) - (w(y,s)) se (_ K 2 >0) \ yt 



< C / |w t | 2 xt, 

where the last inequality follows from (fpQ|). Hence, 

\x ■ (Vw)b x {4} — x • (Vu>) - | 2 xt < C / |t«t| 2 xt. 

To conclude the lemma, we only need to combine (J9)), (jlOl) and (ITT1) . 



(11) 



□ 



3 The projection II and some technical results 

Definition 10. Let II(u,r, X ) be the projection of u into the space of parabolic 
homogeneous caloric polynomials in Q~(X ). In other words 



f 



Q-(X°) 



\D 2 u{X) - D 2 U(u,r,X°) 



hit / \D 2 u{X) - D 2 p{X)\ 2 , 



P^ V2 JQ-(X 



where Vi is the space of caloric polynomials that are homogeneous of degree two 
in space and homogeneous of degree one in time and where we use the matrix 
norm 

i,3 
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Moreover, we will use the notation 

n(n,r,X°) = S(u,r,X°)p UirtX o(X), 
where p u ,r,x° (X) is a second order homogeneous caloric polynomial such that 

\\D 2 p Ui r : x°\\ L °ofQ-} = sup \D 2 p u , r .x«\ = 1, 
Qi 

and S(u,r,X°) £ R+. 

Below are certain properties of the projection that can be easily verihcd. 
Lemma 11. Let u be as in Theorem]]^ Then 

1. II(-,r, X ) is linear, 

2. ||n( M ,r,X°)|| i2(Qn <C\\D 2 u\\ L2{Q - ) forr£ [1,1], 

In i/ie above, C is a constant depending only on the dimension. 

Proof. The first statement follows from the fact that II is a projection. In order 
to prove the second statement we observe that if D 2 u £ L 2 (Q^[) then 

inf / \D 2 u- D 2 P\ 2 < [ \D 2 u\ 2 

and moreover 

\\D 2 u - ^ 2 P|| L2(Qr) > P 2 P|| L2(Qr) - \\D 2 u\\ L2{Q - y 

Thus, 

||5 2 II(u,l,0)|| La(gr) <2|| J D 2 u || i2(Qr) . 
Since _D 2 II(m, 1,0) is constant and II(w, 1,0) homogeneous we can conclude 

||n( u ,i,o)|| L2(Qr) <c\\b 2 u\\ LHQ - y 

The third statement now follows by a simple change of variables. 



□ 



Lemma 12. Assume Hu £ L°°(Q 1 ). Then for every X° — (x°,t°) £ Q l i 2 an d 
r < j, the following inequality holds 



D 2 ( u{rX + X °/ t + t0) -U(u,r,X- } ) 
where C depends only on the dimension. 



L2(Qn <c(h|| il(Q - )+ ||^||^ (Qr) ), 
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Proof. Recall that 



D 2 u 



D 2 u 
d t u 



From Theorem [5] it follows that 

\\D 2 u~ (D 2 u) rtX o\\ L2(Q - {x0)) < C(|M| 11Wr) + ||H M ||^(Qr))r ( " +2)/2 . (12) 
We also observe that 



D 2 u 



M(X°,r) 



-(x°) 
Hu 



Hu 



Hu 



-I 



n + 1 / n + 1 



oi n + 1 



where M(X , r) is a constant matrix with zero trace and / the identity matrix. 
Thus, if 

q X °,r = ~x T [M(X°,r)] nxn x + t[M(X° ' , r)] n+1 , n+1 , 



then 

It follows that 



(D 2 u) r , x o 



b 2 q x o, r = M(X°,r). 

/ b 2 u = b 2 qx ^ r ^ 



Hence, 



L> 2 u - i) 2 g x o : 



< 




'- b 2 u 








Qr(X°) 





< 

L2(Q-{X")) 

Hu 



-I 



n+1 



LHQ-(XO)) 



From this, the definition of II and rescaling the inequality (TT2"j) . the lemma 
follows. □ 



4 The key proposition 



In this section we present Proposition Q3J which is the base of the whole paper. 
First a lemma. 

Lemma 13. Assume that w solves 

Hw(x,t) — f(rx,r 2 t) in Q7 , 

w(x, t) = u r (x,t) — r, 0) on dpQi , 

where 

. . u(rx,r 2 t 
u r (x,t) = — . 

Then 

\\D 2 w\\ L ~ [Q - /2) < C(\\u\\ LHQI) + WDM^-J. 
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Proof. To simplify the notation write S r = S(u, r, 0) and p r — p u , r ,o- Define the 
function 

U = U r — S r p r — (u r — S t P t )q- — X ■ (v (u r — S r p r )Q-j . 

Since Lemma [12l implies 

\\D 2 ( Ur - S r p r )\\ L2{Q - ) < C{ \\ U || £I (Q - ) + || Hu || ioo (g - j ) 

we can use Lemma [9] to obtain 

INI^cq-) < c{M\ L ^) + II^II^(q-))- ( 13 ) 

8 

Now we observe that Hu = Hu so that interior estimates (Lemma |8]) combined 
with (JT3J) imply 

\(u r ~ S r p r ) Q - \ + \x- (V(tl r - S r p r )) Q - | = |Ct r (0)| + |Vu r (Q)| < 

< cGHIlhqd + ll^llL~( Qr ))- 

As a consequence 

||u - S r p r \\ L2{QI) < G(||«|| L i wr) + ||fftt|| LOO(Qr) ). 

To obtain the desired estimate, we use interior estimates (Lemma |5|) on u> to- 
gether with the definition of v. □ 

We recall the notation 

A n Q-| 



X r = 



\Qr\ 

Proposition 14. Let u be a solution to (Qp. Then there exist Cq and C\ 
depending only on the dimension such that if X° e Afl Qi/ 2 an d r < \ then 

gjg^cgn Al/2 > Al/2 

S(r,tt,X")-C 1 (||u|| ilWr) + ||Da«|| i0O(g - ) ) r " 5 ' 

5(r,M,X ) > 2Ci(||u|Ui (Bl ) + \\D 2 v\\ L ~ {Bl) ). 
Proof. For simplicity let X° — and 

u r (X) = w r (X) + S(u, r, Q)p u>rt0 (X) + g r {X), 

where 

Hg r (x,t) = - f '{rx,r 2 t)xK(u(rx,rH)) in Qi , 
g r = ond p Q7, 

and 

Hw r {x,t) — f{rx,r 2 t) in Q7 , 

w r {x,t) = u r (x,t) — S(u,r,0)p urO (X) on d p Q7 . 
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From Lemma [7] and Lemma [T3] it follows that 

L2(Q- ) < C\\f\\ L ™\\xA(u(rx,rH))\\ L 2 ( Q- V (14) 

and 

l|£V[[ £ ~ w - a ) < ^(ll«llzi W r) + H^llL-wr))- ^ 

At this stage we use that D 2 u = dtu = a.e. in {u = 0} and thus 
0= ||£ ,2 Wr|| L 2 (Ai . nQi) = \\D 2 (w r + S(«,r,0)p„, r , + 3r)|| L 2 (ArnQ - ) , 

2 2 

which implies 

\\D 2 S(u,r,0)p Utrfi \\ L2{ArnQ - ) < ||-DV|| i2 (A r nQ7) + IPVllz^nQY)- ( 16 ) 

7 3 7 

From the definition of p u ,r,o it follows that 

||£> 2 2W,o|| L 2 (ArnQI) > ||5 a P»,r,o|| i .( Ar nQ-)V > cA l- ( 17 ) 

2 2 

In addition, the estimates in (fI3)l and ffT5j) imply 

II^Vll i2(Ar .n QI ) < c|l/l[ i -( Ar nor)^ ^ c-II^IL-wr)^- ( 18 ) 

2 

and 

PV|| L2(ArnQT) < C(|| U || Ll(Qr) + ||^|l L » (Ql -))A|. (19) 

2 

Combining {HJ, (HHJ) and (19} we can conclude 

S(u,r,0)Xl < CW&vW^^X? + C(\\u\\ Ll(QI) + \\D 2 v\\l°° {QI ))4> 

from which the desired result follows. □ 

Proposition 15. Let 

Hg = f(x,t)xA inQ~_ k 
\\D 2 9\\ L ^ k) < Ct2-««*W 

where < C 2 and 

\ 2 - } -i < -A 2 -j (20) 
for j = fc,fc + l,fc + 2,...,J. Then there exists a universal constant Co such that 

sup|II(. 9) 2-',0)| <C (C 1 +C 2 ) (21) 
Qi 

for j = k,k+ 1, k + 2, J. 
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Proof. There is no loss of generality to assume that k = 0. In particular, we 
may parabolically rescale g to g(x,t) — g(2~ k x, 2~ 2k t)/2~ 2k . Then g satisfies 
the assumptions with k = and J — k in place of J. If we can prove ([?!]) for 
<7 then it follows for g by scaling back. We may thus assume that k = in the 
proof. 

We may write gj = g\g- , the restriction of g to Q^-j > as the following sum 

3 

ft = I>fc+fc (22) 

fc=0 

where 

Hg = fXA in 

fjj =0 on d p Q~_ 3 

and 

iJ/i fe = in 

hk = fjk-i ond p Q 2 _ k , 

where we, for consistency, identify g_i = g. That gj = 5 in follows by an 
easy induction. It is true, by definition, for j = — 1. If g 3 -\ = g in Q7-j+i then 
i?<?j = iZp in Q2-3 an d 3j = j on dpQ^-j by construction. By the maximum 
principle it follows that gj = g in Q^-j and our induction is complete. 
Next we notice that by ([20]) and Holder's inequality it follows that 



\m j \\ LaiQ -_ j) <c a 2-s\Q-_ i \ 1 / 2 . 

This implies, by Lemma [3 that 

UjWw^iQ-^) < CC 2 2^\Q-_ 3 \ 1/2 (23) 

for some universal constant C. 

Since foj+i is caloric with gj as boundary values it follows from Lemma [7] 
that 

sup |II(^ +1 , 2^,0) | < CC 2 2- j (24) 
QI 

for some universal constant C and all k > j + 1 and all j > 0. 

From (|23[) we may deduce, using Lemma [TT] and a simple rescaling, that 

sup|n(&,2-J',0)|<C7C22-'. (25) 

We need to estimate the projection of ho as well. From parabolic estimates 
(cf. Lemma [7]) we can conclude that 

||£ 2 M L2(Qr) < Cp 2 .g|| i2(Qr) 

and thus, using Lemma [11] again, that 

sup|II(V2- fc ,0)| < CCi (26) 
QI 
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for any k > 0. 

Using the linearity of IT we can finally conclude that for any j < J 

3 

sup \U(g, T\ 0)| < SU P |n(h fc) 2-^,0)1 + \U(~ gj ,2^, 0)| < 
Qi k=o Qi 

<c(c 1 + C 2 ^- k + C 2 )<C (d + 3C 2 ) , 
\ fc=o / 

where we have used (|22[) in the first inequality and (f24|) . (|25|) and ([26]) in the 
second. □ 



5 Proof of the main Result 

The two following lemmata provides us with the result that if S is bounded then 
u is 

Lemma 16 (Quadratic growth implies ). Suppose u is a solution of ([1} 
such that 

sup \u\ < Mr 2 
for all < r < 1/2 and Y° C Q7 n d{u ^ 0}. Then 

2 

\\b 2 u\\ L ~ (Qr) <C(M+ |M| L1(Qn + \\D 2 v\\ L ~ iQI) ) ■ 
Proo/. Take X° <E Q7 and define 

2 

r = sup{r : Q;(X°) n d{u ^ 0} = 0}. 
We now split the proof into two cases: 

Case 1: r < 4. From the definition of r it follows that with 

w = u- v - v(X°) - Vvpf ) • (x° - x), 
then Hw = in Q~(X°). From Lemma[6] 

H^IU(Q r - /2 (x°)) < ^II u IIl~(q,t(xo)) < c (a/ + IP 2 «|| ioo(Q - (x0)) ) , 



and thus 



\D 2 u\\ L ~ {Q;/ a( xo)) < C (M + II^IU (Qn y , 



where C is possibly a larger constant. 
Case 2: r > i. In this case, it is clear that 



w = u- v - v(X°) - Vv(X ) ■ (x° - x), 
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then Hw = in Q k (X°). Lemma [5] then implies 

4 

\\V 2w \\l~{q- { x°)) < c \\ u \\l^q-(xo)) < C I ||«|| £ i (Qr) + \\D 2 v\\ L ^ (Q - [XO)) 
which implies 

\\D 2 u\\l-( Q - /8 (xo)) < C (\H\lhq-) + II^IIl-wd) • 
The combination of the two cases above yields the desired estimate. □ 

Lemma 17 (Bounded S implies quadratic growth). Suppose u is a solution of 
P and let X° e Q7 n d{u ^ 0}. Then for r < 1/4 



sup M < C (s(u,r,X°) + |M| il(Qr) + WHuW^q-, 
>i(x°) 



r 2 . 



The proof is very similar to the proof of Lemma [T3] and therefore we give 
only a sketch of the proof. 

Proof. To avoid cumbersome writing we use the notation 5* = S(u,r, X°) and 

u(rx + x°,r 2 t + t°) 
u r = u r .x° = -jj ' 

where X Q = (x°, t°), throughout the whole proof. The hypotheses of the lemma 
and Lemma [T2l imply 

\\D 2 u r \\ L2{Q - ) < C (s+\\u\\ L1(QI) + \\Hu\\ Loc{Q - } 

Defining 

u = u r - (u r ) Q - - x ■ (v(it r ) Q -J , 
we can reason as in the proof of Lemma [T3] to obtain that 



INL~(Q T ) + IM Q -I + \x ■ (VM) Q - 1 = ||«|| l „ ( q T) + K(o)l + |v« r (o)| 

< c(5+|H| il(Qr) + p^lU (Qr) ) . 

This implies, by the triangle inequality 

IKII £ oc (QT) < C(S + ||«|| £ i C<3r) + WHuW^^-)). 

2 

□ 

We are now ready to give the proof of the main theorem. 
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Proof of Theorem^ In view of Lemma 1171 and Lemma [TBI it is enough to prove 
that 

sup|II(tt,2^X )| < C(\\u\\ Ll{ } + \\D 2 v\\ LOO(Q - } ), (27) 

or 

for j = 2,3,4,... and some universal constant C and every X° £ ^ ^ s 

enough to prove (|2T|) for X° = 0. Once (|27p is proved for X° = a translation 
argument assures that ([27| holds for any X° £ Qy 2 - Then Lemma [T7l provides 
a quadratic bound on the solution which, by Lemma [T5] implies regularity. 
Let us denote by B the set B c N of all j £ N such that 

2-^,0) > 2C ||5 2 «|| £0 o (on + 2d (|H il(Qr) + II^ILco^)) , 

where Co and Ci are as is Proposition [JJ] Naturally if j £ B then 

2-^,0) < 2C p 2 HI^ (Qr) +2d (ll«ll £lWn + II^IU (or) ) 

for all j which in turn implies (|2T|) with C = 2(cq + Ci). Therefore i? consists 
of the "bad" scales where (|2"7)) might not hold. 
If j £ -B then, according to Proposition [TJ] 

-A 2 - 3 > A 2 -j-i . (28) 

Moreover, if {fc, fc + 1, k + 2, J} C B then (g5]) holds for all j £ {fc, k + 
1, J}. In particular, if {fc, fc + 1, k + 2, J} C B and if we split u into 

u = v + g (29) 

where Hv = Hv and i/g = —fxA then 5 satisfies the conditions in Proposition 
[15] for j £ {fc, fc + 1, J}. We have a choice in h and 3 and we may choose g 
such that ||<7|| £ i(Q-) < C|I/IIl°°(q-) by for instance letting g be the convolution 
of the heat kernel and —fxA- Then Theorem [5] implies that 

\\Dg\\ BM o(Q- /2 )<C\\fh^ QI y (30) 

Then it follows, by the triangle inequality, that 

Nlzi(o r )<ll«ll £ i Wl -) + ^ll/llL-(or)- ( 31 ) 

Let j £ N. Then cither j £ B and (|2T|) holds or j £ B and there exists a 
smallest k £ N such that {k,k + l,k + 2, j} £ B. Using ([25]). the linearity of 
II and the triangle inequality we can estimate 

sup |ll(u, 2^,0) -II(w,2- fc ,0)| < sup 2^,0) -n(«,2- fc ,o)| + (32) 
Qi Qi 

+ sup|n( 5 ,2^,o)-n( 3 ,2- fe ,o)|. 

QT 
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Since v S W^ 1 1 (Q 1 ) it follows from Lemma |H] that 

sup |n(s, 2- J ',o)| , sup |n(«, 2-\ 0)1 < c (\\b 2 v\\L°° {Q -) + ll^l^on) (33) 
or QI 

<c(\\D 2 v\\ L ^ QT) + \\u\\ LliQT) ), 

where we used (|3"Tj) in the last inequality. 
We use Proposition [TBI to estimate 

sup \ll(g, 2- -3 , 0) - U(g, 2~ k , 0) | = sup |n(ff, 2~' , 0) | (34) 
Qi QI 

where g = g — IT (5, 2~ k , 0). In particular, by (1501) and Lemma IT21 it follows that 
ll5l!L 2 (g-_j<C||/|| L . (Qr) 2- fc ^ 2 )/ 2 . 

Moreover, by (|2"5|) and our assumption that {fc,fc+l,...,j} C -B the assump- 
tions in Proposition [T5l are satisfied with C\ = C||/Hx,°o(q-) f° r some universal 
constant C. From Proposition [T5l and (fMf it therefore follows that 

sup |n( 5 , 2^,0) - n( 5 , 2- fc , 0)| = sup |n(.g, 0)| < C||/|| Loo(Qn . (35) 
QI Qi 

From (l32l), (1331) and (1351) we can conclude that 



sup |n(«, 2-^,o)| < sup |n(u, 2- fc , o)| + c (\\d 2 v\\l~( Q -) + Mmar)) ■ ( 36 ) 

QI QI 

We need to estimate supg- \H(u, 2~ k , 0)|. Remember that k £ N was the 
smallest constant such that {k, fc + 1, j} C B. This implies that k — 1 ^ B so 

S(u,2- fe +\0) < 2C p 2 «|| ioc(Qr) + 2d (lMI L1(Qn + \\b 2 v\\ L - { Q7)) ■ (37) 
Furthermore, by Lemma [T2l 

||^ 2 (^ - XI(^, 2- fc+1 , 0))|[ i2CQ -_ fc+i) 
< C2 -(„ +2 )( fe -i)/ 2 (||/|| ioo(Qr) + ||K|| il(Qr) ) . 

That is 

sup|n(u-n(M,2- fc+1 ,0),2- fe ,0)| < 
QI 

C]|n (t* - n(«, 2-*- 1 - 1 , 0), 2-*, 0) || Zi - c < 3r) < 

2 -fc-fc(n+ 2 )/2 C || n ( u _ n(M; 2 - fe+1 , 0), 2- fc , 0) || L2(Q - j < (38) 
< 2-( n+2 ) fe / 2 C||L) 2 (w-n(u,2-' £+1 ,0))|| L2(Q - j 



(QjT) ll U HLi(Q-) J , 
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where we used that II(-) is a polynomial in the first inequality, a rescaling in 
the second, Lemma [TT] in the third and Lemma [12] in the last inequality. The 
triangle inequality, ([57)1 and ([551) implies that 

sup|n(u,2- fc ,0)| < 
Qi 



< C2("+ 2 >/ 2 (ll/ll^o-) + IH^q-)) +sup|n(u,2- fc+1 ,0)| < (39) 

<c(\\f\\ L - (QT) + \\u\\ LHQI) ). 

Using dHJ) in flU yields 

sup |n(«, 2-^o)| < c (\\b 2 v\\ L ™ (Q -) + ll«llzi W r)) > ( 40 ) 

Qi 

where we also used that ||/||l°° < C'||£' 2 w||loo . Equation (|40p is valid for all 

3 e N. 

By translating the coordinate system it follows from (|40p that 



sup|n(u,2-^X°)| < C (||D 2 «|| ioo(Qr) + ||u|| il(Qr) ) , (41) 



for any X° G Q^ 2 - 

Lemma [T71 and (|^T|) implies, for all r € (0, 1/4), that 



SUP H < C ( P^H^q-) + ||tt|| il( g-- 

which by Lemma \W\ implies that 

\\D 2 u\\^ (Q - /2) < C (\\D 2 v\\ L ~ iQ ;) + \H\ LHQ -)) , 
and the proof is complete. □ 



6 Regularity of the free boundary 

In this section we prove the second main theorem. The idea is to prove that at 
low energy points (see Definition [TO] below) points and at a scale small enough, 
the solution is non-negative. 

6.1 Weiss' monotonicity formula 

In order to prove the regularity of the free boundary we need to introduce some 
notions from for instance HI. 
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Define the Weiss energy for v(x, t): K' 1 xR — > M to be 
W(r; v, f, (x°,t )) = 1 / (\Vv\ 2 + 2fv+^) G(x, -t)dxdt, 

r JwL"x(-r 2 +t°,t°] V 



t 



and let 

v(rx + x a ,r 2 t + t Q ) 
v r ,(x ,t )(x,t) = ^ ' 

so that 

W(r;v,f, (x ,i°)) = W{l;v n{x o tt o ) ,f r!(x o !t o ) ,0), 

where 

fr,(x°, t) = f(rx + x°,r 2 t + t"). 

Moreover, for a general function u(x,t) define 

Lu = x ■ Vu + 2tu t - 2u. (42) 

The following proposition is a parabolic version of Weiss' monotonicity for- 
mula. Once we know that u enjoys the optimal regularity (Theorem [1]), it can 
be proved in the same manner as in [7] or |12) . 

Proposition 18. Let u be a solution of ([T), ip € C^°(B 3 /4) such that t/> = 1 
on B1/2, and set v = uip. Then there is a continuous function 

F = F(\\f\\ CDi u,\\u\\ L1{Q - ) ,r), 

with F(0) = such that 

W(r;vJ, (x°,t a ))+F(r) 

is a non- decreasing function for < r < 1/2, and in particular for < s < r < 
1/2 there holds 

W(r; v, /, (x°,t )) - W(s; v, /, (x°,t )) + F(r) - F(s) 
\ [ ^^-G{x,-t)dtdxdT. 

T JM«x(-r 2 +t°,t ] ~t 

In view of the proposition above, the limit 

W(0 + ; u, 1) := lim W(r, u, 1) = lim W(r; u, f) 

r— >0 r— >0 

exists. In Lemma 6.2 and Lemma 6.3 in [3], the possible values are determined. 
It turns out that they are in general 0, 15/2 and 15. As is explained in Section 
7.1 in [4] the value 0, corresponding to so-called zero energy points, does not 
occur for X° £ T. 

Definition 19. We say that the the point [x ,t ) G T is a low energy point if 

1 K 

W(0+;uJ, (x°,t )) := lim iW(r;vJ, (x°,t )) = ^. 

From Proposition [18] it follows that the function 
(x,t)^W(0 + ;u,f, (x,t)), 
is upper semi-continuous, and thus the set of regular points is an open set. 
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6.2 The proof 

We can now give the proof of the second main theorem. 

Proof of Theorem^ The proof consists of three steps. The first step amounts 
to prove that if the origin is a low energy point, then u is non-negative close to 
the origin. The second step consists of applying the theory known for the case 
when u has a sign, implying that the free boundary is locally a C 1 graph. These 
two steps prove the first part of the theorem. 

In the third and final step, we observe that if the geometric condition holds 
at the origin, then Lemma 13.3 in [3] implies that the energy is sufficiently low 
(below the threshold 15) in a uniform neighbourhood of the origin. In particular, 
this implies that all points in that neighbourhood are low energy points. Hence, 
the first part of the theorem is applicable in a uniform neighborhood of the 
origin, which implies the second part of the theorem. 

Step 1: For tq small enough, it > in Q~ . We argue by contradiction. If 
this is not true, then there is a solution u of |T]), with 

Ml H q-) + ll/llc™ < M, W(0; u, /, (0, 0)) = y , 

and points (y j ,s 3 ) £ ft 4 (0,0) so that u{y° , s 3 ) < 0. Let (x j ,P) be the 
closest free boundary point, i.e., 

r,j = dist(T, {y\s°)) = dist(V - x\ s j -t j ) = a/R -y j \ 2 + \V - si\. 

Define the rescaled functions 

u(r jX + xj,r?t + V) 

v>(x,t) = • 

j 

Clearly Vj satisfies the equation 

Hvj = f(rjx + x\r]t + t J )x{ Vj7 io}, in Ql_(-a^ , -t J ). 

Moreover, due to Theorem [T] 

sup \vj\ < Cp 2 , for p < -!-, 
Qp lr i 

and by the choice of (y^ , s J ) and (x^ , P), Hvj — f(rjX + x 1 ,r 2 t + P ) in the set 

V - s j ] 

By standard estimates for parabolic equations, we can extract a sub-sequence, 
again labelled Vj , such that Vj — > vq uniformly and 

Hv = X{v Q ^o}i in K n x M~, (0,0) £ T(v ) 



(x,t) : 
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sup |t>o| < Cp 2 , for all p > 0, 
Qp 

and 

Hv a = 1 in {(x,t) : \x - z°\ 2 + \t-T°\ < 1}, v {z°,t°)<Q, (43) 

where 

(z°,t°)= lim (yLl*L ?l^±i) (z°) 2 + r° = l. 
j^-co \ r, r, J 

Observe that ((321) assures that vq ^ 0. Moreover non-degeneracy (Lemma 5.1 
in [¥]) implies that the origin is contained in the free boundary of vq (see also 
section 5.2 in [J]). 

Next we need to use the assumption on the energy functional W(r, u, /, X). 
Since W(r, u, f, X) is uniformly continuous in X for each r > and 

IimW(r,u,/,0) = H, 

the monotonicity formula implies that for each e > there exists an r e > such 
that 

W(r,u,/, (aJ») <^(r, W ,/,0)+e< — + 2e 

if r < r e and j is large enough (J may depend on r). In particular, a rescaling 
implies that 

w(^,v j J{r j x,r 2 j t),Qj<^- + 2e 

if r < r £ is small enough and j large enough. Passing to the limit j — > cx> and 
using that is almost monotone in its first argument we may conclude that 
for any e > 

W(r,Wd,l,0)<y + 2e. (44) 
From (|4~4"]) and the second part of Lemma 9.2 in [3] we can conclude that 

v Q = -{x-e)\, (45) 

for some unit vector e. This is a contradiction to (|4"3"]l since (z°,r°) is, by 
construction, a point in the set Hvq — 1 at unit distance from the free boundary 
so by £45]) v {z°,t°) = | which contradicts (|43|) . 

Step 2: Apply the results from Now we are in the situation of 

Theorem 1.9 in [IT] if we consider u to be defined only in Q~ , i.e., u > in Q~ 
and the origin is a low energy point. Hence, there is a small neighbourhood, 
which might depend on u itself, where the free boundary is a C 1 graph (in the 
parabolic sense). 

Step 3: The geometric condition implies a uniform energy condition. 

Theorem Q] applied to u and then Lemma 13.3 in [4] applied to u(rx,r 2 t)/r 2 im- 
plies that if cr(r) is large enough, r small enough (both depending on IMI^ijq-) 
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and ||/|| C Din i(Qr) ) and 

MD ({x : u{x, -r 2 ) = 0} n B r ) 

> <?{r), 

r 

then 

W(p;u,f,X) < 15 -s , 

for all X G Q~i 2 , p < r small enough (depending on j and ||/ ||c*Dini(Q-))- 

In particular, all free boundary points in Q~j 2 are ^ ow ener gy points. Hence, we 
can apply Step 1 and Step 2 to conclude that in Q~, 2 , the free boundary is a 
C 1 graph. 
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